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A New Differential Absorption Lidar to Measure Sub- 
hourly Fluctuation of Tropospheric Ozone Profiles in the 
Baltimore - Washington D.C. Region 

Abstract. Tropospheric ozone profiles have been retrieved 
from the new ground based National Aeronautics and 35 
Space Administration (NASA) Goddard Space Flight Center 
TROPospheric OZone Differential Absorption Lidar (GSFC 
TROPOZ DIAL) in Greenbelt, MD (38.99°N, 76.84°W, 57 
meters ASL) from 400 m to 12 km AGL. Current atmo- 
spheric satellite instruments cannot peer through the op- 
tically thick stratospheric ozone layer to remotely sense 
boundary layer tropospheric ozone. In order to monitor this 
lower ozone more effectively, the Tropospheric Ozone Li- 40 
dar Network (TOLNet) has been developed, which currently 
consists of five stations across the US. The GSFC TROPOZ 
DIAL is based on the Differential Absorption Lidar (DIAL) 
technique, which currently detects two wavelengths, 289 
and 299 nm. Ozone is absorbed more strongly at 289 nm 
than at 299 nm. The DIAL technique exploits this differ- 45 
ence between the returned backscatter signals to obtain the 
ozone number density as a function of altitude. The trans- 
mitted wavelengths are generated by focusing the output of a 
quadrupled Nd:YAG laser beam (266 nm) into a pair of Ra- 
man cells, filled with high pressure hydrogen and deuterium. 50 
Stimulated Raman Scattering (SRS) within the focus gen- 
erates a significant fraction of the pump energy at the first 
Stokes shift. With the knowledge of the ozone absorption co- 
efficient at these two wavelengths, the range resolved number 
density can be derived. An interesting atmospheric case study 55 
involving the Stratospheric-Tropospheric Exchange (STE) of 


ozone is shown to emphasize the regional importance of this 
instrument as well as assessing the validation and calibration 
of data. The retrieval yields an uncertainty of 16-19% from 
0-1.5 km, 10 - 18% from 1.5 -3 km, and 1 1 - 25% from 3 km 
to 12 km. There are currently surface ozone measurements 
hourly and ozonesonde launches occasionally, but this sys- 
tem will be the first to make routine tropospheric ozone pro- 
file measurements in the Baltimore-Washington D.C. area. 


1 Introduction 

Tropospheric ozone above the ground level has been histor- 
ically difficult to measure directly due to its relatively short 
lifetime and nonlinear formation (Stevenson et al. (2013)). 
It is an important greenhouse gas, pollutant, and source of 
OH radicals. Its contribution to global warming from the 
preindustrial era to the present is regarded as the third most 
important, following those of carbon dioxide (CO 2 ) and 
methane (CH 4 ) (IPCC (2007)). Ozone is also toxic to hu- 
mans and vegetation because it begins oxidizing biological 
tissue and causes harmful respiratory effects in instances 
of long exposure (McDonnell et al. (1999)). Tropospheric 
ozone within the Planetary Boundary Layer (PBL) in the 
Baltimore-Washington D.C. area is typically dominated by 
photochemistry with local ozone precursor emissions, but 
within the free troposphere there are other source classifica- 
tions such as 1) Stratospheric-Tropospheric Exchange (STE), 
2) advection from upwind sources and 3) regional redistri- 
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bution from convection and/or lightning (Thompson et al. 
(2008)). 

During an STE event, large concentrations of ozone are 115 
exchanged between the stratosphere and the troposphere dur- 
ing certain meteorological events, such as an unusually low 
Polar or Subtropical jet stream (Langford (1999)). A recent 
analysis of summertime tropospheric ozone profiles in the 
Baltimore-Washington D.C. area showed that STE accounted 120 
for greater than 50% of the free tropospheric ozone column 
on 17% of days sampled (Yorks et al. (2009)). The second 
source of ozone in the free troposphere occurs when ozone 
is transported aloft or advected with varying layer thickness 
and concentration downwind of the major ozone precursor^ 
production sites, potentially resulting in an ozone exceedance 
for rural and less populated areas (NCADAC (2013)). The 
subsidence of the ozone layer aloft has been shown to cause 
air quality issues for a region that may not have actually pro- 
duced the ozone or its precursors (Langford et al. (2010)). 130 
The third source of ozone in the free troposphere is not a pri- 
mary source but rather the ozone precursor NOx is induced 
by lightning flashes, which accounts for 5 - 10% of the over- 
all tropospheric ozone budget (Schumann and Huntrieser 
(2007), Yuan et al. (2012)). 135 

With such dynamic sources of ozone in the troposphere, 
it is essential to monitor the vertical distribution of tropo- 
spheric ozone. Although the vertical location of atmospheric 
trace gases (such as CO 2 ) are remotely monitored by satel- 
lites, these measurements in the ultra-violet (UV) regime, uo 
especially near the surface, are difficult to obtain due to 
the optically thick stratospheric ozone layer strongly atten- 
uating the signal (Fishman et al. (1990)). Ozone profiles 
from the Tropospheric Emission Spectrometer (TES) and 
the Ozone Mapping Instrument (OMI) have been reported, us 
but have usually required a large amount of prior knowl- 
edge and/or ancillary satellites (such as NASA’s Microwave 
Limb Sounder (MLS)) are needed to determine the residual 
amount of ozone in the troposphere (Ziemke et al. (2006), 
Bowman et al. (2002)). Even if the tropospheric ozone profile 150 
is reported, it may be a seasonal or monthly mean, in which 
the variability of ozone gives the profiles marginal statistical 
value (Thompson et al. (2003)). For ozone profiles that are 
reported for a direct overpass, due to the large footprint of the 
satellite, this method may not be able to resolve small scale 155 
fluctuations of tropospheric ozone concentration, especially 
during quickly evolving pollution episodes. By monitoring 
and analyzing tropospheric ozone from a ground based in- 
strument, such as an ozone lidar, the signal is not attenuated 
via stratospheric ozone and the tropospheric ozone concen - 160 
trations can be detected at a much higher temporal and verti- 
cal resolution. 

For these reasons, the National Aeronautics and Space 
Administration (NASA) Goddard Space Flight Center TRO- 
Pospheric OZone Differential Absorption Lidar (GSFC i6s 
TROPOZ DIAL) has been developed to take routine 
measurements of tropospheric ozone near the Balti- 


more - Washington D.C. area (Greenbelt, MD 38.99°N, 
76.84°W, 57 meters ASL). This instrument has been 
developed as part of the ground-based Tropospheric 
Ozone Lidar NETwork (TOLNet), which currently con- 
sists of five stations across the United States (http://www- 
air.larc. nasa.gov/missions/TOLNet/). The primary purpose 
of the instruments within TOLNet is to provide regular, 
high fidelity profile measurements of ozone within the tropo- 
sphere, particularly the changing ozone dynamics and lam- 
inae inside the PBL. The network will also serve to vali- 
date the upcoming GEOstationary Coastal and Air Pollution 
Events (GEO-CAPE) mission, scheduled to be launched in 
the 2020 timeframe. Another objective of TOLNet is to iden- 
tify a brassboard ozone lidar instrument that would be suit- 
able to populate a network to address the needs of NASA, 
the Environmental Protection Agency (EPA), and the Na- 
tional Oceanic and Atmospheric Administration (NOAA) 
air-quality scientists and managers who increasingly express 
a desire for ozone profiles. 

Similar ground based instruments have been developed 
in TOLNet (McDermid et al. (2002), Alvarez et al. (2011), 
Kuang et al. (2013)), but the GSFC TROPOZ DIAL is the 
first to measure ozone profiles in the Mid-Atlantic region. 
Aside from the site at the Jet Propulsion Laboratory’s Ta- 
ble Mountain Facility (JPL TMF), this is the only instrument 
currently utilizing high-pressure Stimulated Raman Scatter- 
ing (SRS) cells for wavelength generation (Haner and Mc- 
Dermid (1990)). Because the remote JPL TMF site is in the 
Angeles National Forest at an altitude of 2.3 km ASL, the 
GSFC TROPOZ DIAL will be the first system to use this 
technique to determine quantitative ozone profiles near sea 
level directly in an urban environment. 

At the surface, in situ UV analyzers and photometers are 
the most common form of ozone monitors. These instru- 
ments utilize photometry and the known extinction coeffi- 
cient at a specific UV wavelength (usually 254 nm with a 
Mercury vapor lamp) and the Beer-Lambert Law to derive 
the known ozone concentration in a volume of sampled gas. 
The reduction in light intensity between the two photodetec- 
tors over a known sampled path length is directly related to 
the ozone concentration. Most monitors require calibration 
with a zero air source and are usually accurate to within 5% 
of the known value (Puzak (1987)) and can report values with 
very high temporal resolution. 

The EPA has set up a network of surface ozone monitoring 
and detection that covers a large spatial and temporal region 
in the last few decades (Demerjian (2000)). These ozone con- 
centrations are typically reported in hourly average values or 
8 hour average values and are useful for analyzing surface 
trends in ozone. As mentioned earlier ozone can be trans- 
ported aloft and advected, preventing the surface monitors 
from truly representing the ozone distribution in that region. 
The vertical distribution of ozone is an important parameter 
in analyzing daily and seasonal trends and events that sur- 
face monitors cannot provide. If there is a dramatic change 
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in ozone at the surface that cannot be explained with the in 
situ data, it is essential to analyze the vertical profile of ozone 
above the surface as a function of time.. 225 

The common historical method for vertical ozone detec- 
tion is sending a balloon-borne electrochemical concentra- 
tion cell (ECC) instrument through the atmosphere, but this 
may not be on the continuous scale that is necessary to fully 
characterize small scale fluctuations in ozone. The ECC pro- 230 
files ozone with a 100 m vertical resolution from the sur- 
face to 35-km altitude with an accuracy of 10% for tropo- 
spheric mixing ratios, except in the case of very low mix- 
ing ratios (<10 ppb) when the accuracy may be degraded to 
15% (Newchurch et al. (2003), Komhyr et al. (1995), Olt -235 
mans et al. (1996)). 

Due to complex formation and transport of ozone it is typ- 
ically necessary to measure for extended periods of time to 
fully characterize an ozone pollution episode, which favors a 
continuous measurement as opposed to an isolated or several 240 
profiles. Using lidar techniques in mostly cloud free condi- 
tions, it is possible to measure the continuous vertical dis- 
tribution of the dynamic atmosphere by detecting the return 
signal as a function of altitude. A stationary, continuous time 
series of vertical profile measurements can be obtained more 245 
cost effectively with a DIAL instrument such as described in 
this paper than with the launch of multiple ozonesondes. 

There have been successful DIAL instruments that have 
utilized the ground based monitoring method, dating back 
to the 1970’s, with a focus on water vapor (Browell et al. 250 
(1979)) and temperature profiles (Murray et al. (1980)). 
These concepts were then successfully modified and demon- 
strated to yield accurate results for ozone detection (Uchino 
et al. (1978), Megie et al. (1985), Browell (1989)), with 
specific emphasis on monitoring and maintaing a historic 
record of ozone trends in the stratosphere (McDermid et al. 
(1991), McGee et al. (1993), McGee et al. (1995)). Be- 255 
cause tropospheric ozone can create air pollution episodes, 
it was necessary to develop DIAL instruments, such as the 
GSFC TROPOZ DIAL, and a network, such as TOLNet, to 
make routine measurements of the ozone concentrations in 
the lower free troposphere and the PBL, which can fluctuate 
much more rapidly that that of the stratospheric ozone. 260 

Similar instruments have been deployed as nadir viewing 
airborne systems in previous NASA (Browell et al. (1983), 
Browell et al. (1987), Browell et al. (1993)) and NO A A 
(Senff et al. (2010)) missions. Although there have been suc- 
cessful measurements of tropospheric ozone profiles using 265 
these instruments, there is a sizable operational cost increase 
and logistical preparation for aircraft systems over ground 
based instruments (Ancellet and Ravetta (1998)). Although 
an airborne experiment has many advantages, especially with 
the use of onboard in situ monitors and remote sensing instru- 270 
ments, it can be difficult to quantify ozone near the surface 
due to aviation restrictions. Furthermore, past efforts to pre- 
dict and monitor high levels of ozone in large cities (Houston 
(Banta et al. (2005)) and Los Angeles (Menzies and Shumate 


(1978)) from airborne missions have been hampered by the 
lack of information of the ozone profile and the extent of 
vertical mixing in the first few kilometers above the surface 
(Morris et al. (2010)). 

The metropolitan area surrounding the GSFC TROPOZ 
DIAL is affected by locally produced pollutants from highly 
populated urban areas as well as ozone precursors advected 
from the Ohio River Valley/Appalachian mountain regions. 
Depending on the meteorological conditions, these pollu- 
tants can cause non-localized unfavorable air quality issues 
for smaller rural areas as well as the Chesapeake Bay water- 
shed to the south and east. If the eastward propagation of pol- 
lutants is inhibited, a stagnant air mass over the Baltimore- 
Washington D.C. area results, and is the ideal condition for 
a high pollution episode (Delgado et al. (2014)). During 
these forecasted high ozone episodes, the GSFC TROPOZ 
DIAL has been designed to move from routine monitoring 
to continuous monitoring. The Mid-Atlantic region displays 
a moderate amount of seasonal ozone variability with the 
highest concentrations occurring in the late July time pe- 
riod (Lehman et al. (2004)). Because of this, during the fall 
and winter seasons, the TROPOZ focuses its observations on 
monitoring ozone in the free troposphere and the aforemen- 
tioned source contributions, such as STE. This paper investi- 
gates the steps necessary to construct the GSFC TROPOZ 
DIAL in order to retrieve and better characterize the ver- 
tical distribution of tropospheric ozone in the Baltimore- 
Washington D.C. area, in which a regionally important vali- 
dation test case of an unusually low tropopause, resulting in 
the exchange of high amounts of ozone between the strato- 
sphere and the lower free troposphere, or STE, is presented 
and analyzed. 

2 The Differential Absorption Lidar (DIAL) Method 

Fundamentally, the DIAL method relies on the fact that the 
ozone molecule has a strongly varying spectral absorption 
curve in ozone’s Hartley band region. Figure 1, adapted 
from Malicet et al. (1995), illustrates the ozone absorption 
cross section ( ao 3 ) and the TROPOZ DIAL wavelengths 
that were chosen for the large gradient in ozone absorp- 
tion (A<to 3 ), which will yield usable atmospheric informa- 
tion. Different colors in the figure correspond to the tem- 
perature dependence of the absorption cross sections. This 
method is optimized when the ratio of the difference in the 
ozone absorption cross section to the difference in wave- 
length ( A<to 3 / AA) is maximized within certain atmospheric 
constraints. 

Retrieval of useful ozone profiles during daytime hours re- 
quires that the instrument utilizes wavelengths less than 300 
nm to avoid saturation of the detectors from sunlight, and 
to maximize the signal to noise ratio (SNR) of the system 
(Proffitt and Langford (1997)). Although the ozone absorp- 
tion curve continues to rise sharply as the wavelength de- 


275 

280 

285 

290 

295 

300 

305 

310 

315 

320 


4 


Sullivan: Tropospheric Ozone DIAL in the Baltimore-Washington D.C. Region 


creases (increasing the A<to 3 /AA), the signal is attenuated 
much too strongly to yield usable information above above 
the PBL. Essentially, shorter wavelengths are absorbed so 325 
strongly by ozone that the return signal would be unusable 
above the PBL, especially during pollution episodes. These 
issues constrain the choice of wavelengths to between 280 
and 300 nm. This ensures the system has a large enough sig- 
nal to easily retrieve ozone profiles through the PBL and into 
the lower free troposphere. 

Lrom Ligure 1, 288.9 nm has an absorption cross section 
nearly 3 times as large as 299.1 nm and both wavelengths 330 
can be generated reliably (and much more cost-effectively 
than dye or tunable cavity lasers) using SRS in high pres- 
sure Raman cells (Haner and McDermid (1990)). The trans- 
mitted wavelengths are generated by focusing the output of 
two frequency-quadrupled Neodymium-doped Yttrium Alu- 
minum Garnet (Nd:YAG) laser beams (266 nm), each into 335 
a pair of high pressure Raman cells, one filled with hydro- 
gen and the other deuterium. SRS within the focus generates 
a significant fraction of the pump energy at the first Stokes 
shift. With the knowledge of the ozone absorption coefficient 
at these two wavelengths, the range resolved ozone number 340 
density can then be derived. The more absorbing wavelength, 
at 288.9 nm, is called the “on” wavelength and the less ab- 
sorbing wavelength, at 299.1 nm, is called the “off' wave- 
length. This process of SRS and the generation of the nec- 
essary DIAL wavelengths will be discussed more in the Ex- 
perimental Constmction section of the paper. Because the in- 
elastic Raman conversion is done prior to transmission into 
the atmosphere, and detection of atmospheric scattering is 
done at these same transmitted wavelengths, the TROPOZ 
retrieval used is an elastic technique. 345 

2.1 Introduction to Elastic Lidar 


P° n ( r ) _ Pon{ r ) [-2 fg(ao n (r)-a of f(r))dr] (O') 

Poff(r ) (l off (r) 

Assuming ozone has negligible scattering, the extinction co- 
efficient can be regarded as the absorption coefficient. The 
absorption coefficients, assuming that the sampled volumes 
have the same number density of ozone, can be written as, 

^o//(^) ^o 3 A(Jq 3 ^^-atmospherei (3) 

where N 03 is the ozone number density and Aero 3 is the 
difference in corresponding ozone absorption cross sec- 
tions taken for the two wavelengths in Ligure 1. The term, 
A otatmosphere, is the difference in the scattering properties 
of the atmosphere between the two DIAL wavelengths. This 
includes spectral differences in Rayleigh scattering, aerosols, 
and other interfering gases that have similar absorption prop- 
erties in the “on” and “off’ wavelengths. 

Linally the DIAL equation can be written in terms of the 
ozone number density as, 

at ( \ ^ ^ r r Pon(f) /5o//(t"), A CXatmosph ere 

° Ar) ~2Kao 3 Tr n[ P of f(r) /3 on (r) l+ Aoo7~ 

(4) 

In reality, lidar return signals are not recorded or analyzed as 
continuous functions, but rather as values in discrete range 
bins, A r. It is then possible to write the discrete DIAL equa- 
tion (Milton and Woods (1987)) in terms of the range bins 
specified, other interfering gases (IG) and aerosol corrections 
as, 


The general lidar equation, taken from Middleton and 
Spilhaus (1953), is for elastic backscatter in a homogeneous 
volume of the atmosphere. The elastic lidar equation is, 


N, 


o 3 


1 ln( - p °//( r ' + Ar ) AnQ) 

2Aero 3 Ar Poff{r) ’ P on (r + Ar) 


C)-D 

(5) 


P(r) = P( °)^f^)e l “ 2/ » a(r>r ' 1 + A, (1) 

where P(r) is the power received from range r, P( 0) is the 350 
average transmitted power during the laser pulse at the sur- 
face, 77 is the receiver efficiency, A is the receiver area, r is 
the range to the scattering volume, c is the speed of light, r 
is the laser pulse duration, /? and a are, respectively, the at- 
mospheric backscatter coefficient and atmospheric extinction 
coefficient at range r, and Pj, is the background radiation at 
the transmitted wavelength. 

Writing equation (1) for P on (the lower and more absorbed 355 
wavelength) and P D // (the higher and less absorbed wave- 
length) and taking their ratio, many of the system constants 
drop out. Assuming the solar background can be quantified 
and subtracted, this can be rewritten as 


where, 

fe//(r + Ar) 0 m (r) 

Poff{r) /3 on (r + Ar) 

and 

_ A a m oi Aa aer NiqAcyig 

“ ~Acto7 + A^7 + Aao 3 ' () 

Equation (5) is referred to as the DIAL equation and it is of 
great interest because it is a self calibrating technique that can 
determine the number density of ozone with only the known 
ozone absorption cross sections and the power returned to 
each detector. 

The correction term C is due to the derivative in the 
amount of photons backscattered to the detector from the 
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ambient atmosphere. This can be a combination of Rayleigh 
scattering and aerosol scattering if aerosols are present. 
Rayleigh extinction (which is the first term in D) is typically 
based on a A -4 scaling, which develops a slight difference 
in the molecular extinction that is fairly easy to interpolate 415 
(Measures (1983)) with additional meteorological terms for 
corrections given by a reference standard atmosphere (U.S. 
Standard (1976)). With the knowledge of the Rayleigh ex- 
tinction values, the Rayleigh backscatter term in C is com- 
puted using the assumed Rayleigh phase function. 420 

The second term in D is due to the presence of aerosols 
and the wavelength dependence of the aerosol extinction co- 
efficient. Depending on the optical properties of the aerosols 
in the atmosphere during the measurement, this value can 
change as well, leading to an additional necessary correc -425 
tion (Bosenberg (1998)). In most cases the correction for 
the additional backscattered light from aerosol gradients to 
the detector is much larger than the necessary correction 
for the aerosol extinction correction at these wavelengths 
(Weitkamp (2006)). The aerosol connection term in D is de -430 
termined using an iterative method, in which using various 
valid assumptions it is possible to determine the necessary 
ozone number density correction. The correction is then cou- 
pled with an assumed lidar ratio to determine the correction 
factor due to aerosols in C. This will be discussed further in 435 
the Data Processing section of this paper. 

The third term in D is due to the interference of external in- 
terfering gases (IG), such as sulfur dioxide (SO 2 ) and nitric 
oxide (WO 2 ). In the wavelength region used for the present 
work, (SO 2 ) has relatively similar absorption cross sections 440 
to ozone, but the number density, (Nig) = ( Nso 2 ), is typi- 
cally so low that the contribution it makes to the calculation 
of D is small compared to that from the Rayleigh or aerosol 
terms. Nitric oxide (NO 2 ) has similar extinction cross sec- 
tions which increase with respect to wavelength, but for the 
current DIAL wavelengths, the result should be reasonably 
unaffected (Kovalev and Eichinger (2004)). 


SRS is a direct two-photon, third-order nonlinear opti- 
cal process, which results in inelastic light scattering with- 
out the requirement of phase matching (Mollenauer et al. 
(1992)). This process is efficient for gases such as hydrogen 
and deuterium because they have high-frequency vibrational 
modes and narrow spectral linewidths (Shen and Bloember- 
gen (1965)). Multiple orders of Stokes frequencies are sub- 
sequently generated, but for use in this experiment, only the 
first Stokes shift in each Raman gain medium is useful. 

In order to maximize the generation of the first Stokes fre- 
quency, a low numerical aperture pumping geometry and a 
long optical interaction length (Carlsten et al. (1984)) has 
been established which is the reason for the length of the 
Raman cells (76 cm) and the focal points (f = 1 m) of the 
lenses. Because this process is a photon to photon conversion 
method it is also favorable to have the largest volume avail- 
able to increase efficiency of the SRS (Falsini et al. (1985)). 
A low-angle pumping geometry can also reduce the angle 
condition that is auspicious for phase matched four wave 
mixing (Duncan et al. (1986)). 

This difference in frequency emitted and pumped is di- 
rectly related to the energy difference between the first ex- 
cited state and the ground state in the Raman active gas, 
which is known as the Raman frequency and is adapted from 
Measures (1983) in Table 1. Therefore, the Raman frequency 
of the Raman active gas is independent of the absolute energy 
of the pump photon and the energy of the generated photon 
is a function of the pump frequency and the species used as 
the Raman active gas. In the steady state, this process yields 
the continuous generation of coherent light pulses with the 
same phase as the pump photons. 

The steady-state gain coefficient for the vibrational stim- 
ulated Raman scattering is related to the photon conversion 
efficiency in the system and can be written as 


2Nj\ s 3 da 
nhcT duj ’ 


( 8 ) 


3 Experimental Construction 

3.1 Raman Cells and Stimulated Raman Scattering 

(SRS) 

450 

Although previous work has shown accurate tropospheric 
ozone retrievals with dye laser sources (Kuang et al. (2011)), 
and with Ce:LiCAF in a tunable cavity (Senff et al. (2010)), 
they are expensive and typically delicate systems. The lasers 
are also required to undergo frequent maintenance proce -455 
dures and rigorous optimizations during normal operation 
of the crystals and cavities due to small temperature fluc- 
tuations. For these reasons it was advantageous to move to- 
ward the use of a stable Nd:YAG laser with the known fourth 
harmonic frequency pumping a pair of high-pressure Raman 460 
cells to generate the transmitted frequencies. 


where X s is the first Stokes wavelength, Nj is the Boltzmann 
number density, V is the Raman full line width, ^ is the dif- 
ferential Raman scattering cross section, h is Planck’s con- 
stant, and c is the speed of light (Komine (1986)). 

For this to be an efficient technique, the steady-state gain 
must be at an appropriate threshold level to induce the first 
order Stokes frequencies, while suppressing higher order fre- 
quencies. A straightforward method to increase the number 
density, Nj, is to increase the pressure of the gas in the 
confined Raman cell volume. Although this does increase 
the steady state gain, it has an adverse effect on the over- 
all conversion efficiency for the first Stokes shift because the 
gain grows so large that it begins to activate other higher or- 
der processes, collision effects, motional narrowing and four 
wave matching (Witkowicz and May (1976), Dion and May 
(1973)). 
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One other important factor is that the Raman full line 515 
width, r, is linearly dependent on the total pressure of the cell 
in the pressure region of interest (Murray and Javan (1972)). 
Because of this, the steady state gain and photon conversion 
efficiency are more efficient at higher pressures due to the 
broadening Raman full line width. 520 

Therefore, to obtain the benefits of the high-pressure 
regime, while minimizing the effects of an unnecessarily 
high Raman active number density, a buffer gas has been 
added to the Raman active gas. Helium has been used, as 
it is not a Raman active medium in the wavelength regime of 525 
ozone detection. 

3.2 Raman Cell Characterization 

There is a discussion of this method of using a Nd:YAGs3o 
laser and Raman cells to produce appropriate wavelengths 
for ozone detection in Haner and McDermid (1990), which 
was the starting point for the empirical Raman tests that were 
performed to optimize each Raman cell. Because of the sen- 
sitivity of the photon conversion efficiency in laser beam 535 
characteristics, such as pump beam diameter and confocal 
beam parameter entering each Raman cell, it was necessary 
to test the cells for various pressures and buffer gas concen- 
trations. Unlike other instruments described in the literature, 
this instrument is housed inside a portable trailer, which has 540 
limited the optical interaction length of the Raman cells to 
1.8 m. Due to these constraints, the Raman cells were com- 
pletely designed and fabricated internally for the specific ap- 
plication for use in the GSFC TROPOZ DIAL. 

The simplest way to analyze the photon efficiency was to 
adapt a dispersive pellin broca prism to separate the wave- 
lengths into discrete spots and subsequently view the beam 545 
with a fluorescent card. This provided the ability to quali- 
tatively determine which pressure combinations suppressed 
higher order Stokes shifts and which magnified higher or- 
ders. Power measurements of the beams produced by the var- 
ious Stokes shifts then allowed a quantitative selection of the 550 
optimal partial pressures for each of the Raman active gases. 
Because this is a direct two-photon conversion process from 
the pump to the respective first Stokes, it was necessary to ob- 
tain a baseline for photon conversion by measuring the power 
of the pump photons, at 266 nm, passing through an evacu- 555 
ated cell. 

The empirical results for hydrogen and deuterium are 
shown in Figure 2a and Figure 2b. These plots are shown 
as photon conversion efficiency (left y-axis), residual pump 
transmission (right y-axis) and the total cell pressure (x-axis). 560 
Due to safety concerns regarding the design of the Raman 
cells, the maximum pressure tested was 1000 psig. Multi- 
ple power measurements were taken and each data point cor- 
responds to the mean with uncertainty bars of one standard 
deviation. Initially, power measurements were made without 565 
a buffer gas, but the photon conversion efficiencies were so 
low ( 10 % or less) that they are not reported in this figure. 


Once the buffer gas was introduced into the Raman cell, the 
efficiencies began to improve. The three colors in Figure 2a 
and 2 b correspond to three different pressure combinations 
in each cell. The line style corresponds to the wavelength 
of interest that was measured, which was either the pump 
(PUMP), the first Stokes (SI), or the second Stokes (S2). 

For example in Figure 2a, the red curves correspond to a 
partial pressure of 200 psig of hydrogen and is varied with 
an increasing presence of helium to 900 psig. The solid line 
is the pump wavelength, which steadily grew in the experi- 
ment as the partial pressure increased. The dotted line is the 
SI which grew steadily, but dropped after 700 psig total cell 
pressure. The dashed line S2, slowly decreased throughout 
the experiment. From these results it can be inferred that with 
an increasing amount of buffer gas will suppress the the sec- 
ond Stokes wavelength and increase the output power of the 
first Stokes and residual pump wavelengths. 

After these tests were completed, it was determined that 
a pressure combination of 200 psig H 2 /600 psig He would 
be the most efficient (53%). Because 266 nm will hope- 
fully be used for future work, it was important to maintain 
some residual amount of the transmitted pump beam. A sim- 
ilar process was completed for the deuterium cell, which the 
pressure combination of 300 psig D 2 /50Q psig He would be 
the most efficient (27%). The difference in these two efficien- 
cies can be largely attributed to the fact that the hydrogen Ra- 
man scattering cross section, (^) from ( 8 ), is almost twice 
as large as the deuterium scattering cross section. 

3.3 System Description 

The schematic of the GSFC TROPOZ DIAL system with the 
hydrogen (H 2 ) and deuterium (/J 2 ) Raman cells installed in 
a portable 13 m trailer is depicted in Figure 3 and an im- 
age of the trailer is in Figure 4 below. Lidar systems have 
been shown to produce comparable results in mobile units as 
compared to traditional stationary laboratories, but have the 
advantage of being transported if an event of interest arises 
(McGee et al. (1991)). 

A Spectra Physics Quanta-Ray pulsed Nd:YAG laser has 
been optimized with two independent parallel laser cavities. 
These are each comprised of an oscillator/preamplifier with 
10 mm YAG rods, followed by two power amplifiers. The 
beams then pass through potassium dideuterium phosphate 
(KD*P) crystals. This non-linear optic medium serves the 
purpose of generating the second and fourth harmonics of 
the fundamental Nd:YAG frequency. The laser has been op- 
timized for the conversion of the fundamental to the fourth 
harmonic at a wavelength of 266 nm. These initial pump 
beams are each approximately 9.6 mm in diameter with a 
beam divergence of 0.42 mrad. 

The two emitted beams are each steered to the Raman cells 
using two turning mirrors which have been coated to opti- 
mize reflectance at 266 nm. The beams pass through con- 
verging lenses (LI, L2), with focal length 1.0 m, which fo- 
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cus each beam waist near the center of the 1.8 m long Raman 
cells. By focusing the beam waist near the center of the Ra- 
man cell, the optical interaction volume is maximized which 625 
will give the highest possible conversion of pump photons 
into 1st Stokes shift photons (McDermid et al. (2002)). 

Both ends of the Raman cells are sealed with sapphire opti- 
cal windows 1.2 cm thick and 3.8 cm in diameter; the flanges 
used to hold these optics allow for a final optical aperture of 630 
2.2 cm. The sapphire was chosen in preference over fused 
silica as the window optical material because it is more ther- 
mally and mechanically durable under conditions of contin- 
uous laser illumination. After the beams exit the Raman cell 
they are collimated with a 1.0 m re-collimating lens (L3,L4) 
to minimize beam divergence. The beams are then steered 
and transmitted through the hatch (90 cm x 92 cm) of the 635 
trailer and into the atmosphere. 

The GSFC TROPOZ DIAL has a large 45 cm diameter 
Newtonian telescope for detecting 289/299 nm at higher alti- 
tudes (High) in the free troposphere and four smaller 2.5 cm 
telescopes to obtain signal near the surface (Middle, Low) for 
289/299 nm (McGee et al. (2005)). For the 45 cm telescope, 640 
the lidar return is focused at the field stop with a 1.0 mrad 
aperture, which defines the field of view (FOV) of the tele- 
scope. The 2.5 cm telescopes have a much wider FOV, 4.5 
mrad, which can accommodate the large dynamic range of 
the Rayleigh backscattered laser light in the near field. There s« 
are very narrow band (<1.0 nm FWHM) interference filters 
(IF 289,IF 299) in front of all of the detectors to enhance the 
amount of usable signal and decrease the amount of ambient 
background light. These filters are critical to the system be- 
cause they are spectrally narrow enough for signal detection eso 
during the daytime, when tropospheric ozone is most active. 
After the backscattered light passes through the interference 
filters it is detected by photomultiplier tubes (PMTs). Neu- 
tral density filters have been placed to avoid saturation of the 
PMTs. Electronic gates are also used on the PMTs to avoid 655 
saturation. 

Individual photons are counted by a transient recorder (LI- 
CEL TR 20 - 80, GmbH, Germany) at a maximum counting 
rate of 300 MHz. This process is controlled by Lab View soft- 
ware which was originally developed for use in other NASA 660 
lidar instruments and modified to fit the needs of the GSFC 
TROPOZ DIAL. The system is triggered via delay genera- 
tors to interleave the two laser pulses (10 ms offset) to reduce 
crosstalk between the two channels, which has been shown to 
decrease electrical system noise in the signal by shifting theses 
timing of this noise beyond the temporal region of data ac- 
quisition. All of the important specifications for hardware for 
each altitude and wavelength channel are listed in Table 2. 

The large 45 cm telescope is mounted so that its axis of 
direction is fixed and can be re-directed only by raising and f .70 
lowering the hydraulic legs of the trailer. The small 2.5 cm 
telescopes, however, are set in adjustable mounts so that they 
can each be directed independently. The alignment of the sys- 
tem is accomplished by steering the outgoing beams into the 


region of the 45 cm telescope’s FOV and optimizing the re- 
turn signals in the channels associated with this telescope. 
Optimizing these signals entails maximizing the signal level 
in the far field. After this alignment is achieved, the axial di- 
rections of each of the four smaller telescopes are adjusted 
so that the return signals associated with them are optimized 
in the same manner as described for the large telescope. This 
process has been done during cloud free nighttime hours to 
minimize solar background and potential contamination of 
clouds and interfering gases. Geometrically, the 45 cm tele- 
scope should be fully overlapped at 462 m AGL and the 2.5 
cm telescopes should be fully overlapped at 105 m AGL. 

4 Data Processing and Algorithm 

The data acquisition system in the GSFC TROPOZ DIAL 
collects 20 second (1000 laser shots) integrations of the de- 
tected backscattered light. This is a low enough temporal 
resolution to make the assumption of a homogeneous atmo- 
sphere during the time of each individual profile. These 20 
second integrated profiles will be referred to as the raw signal 
in this paper. To obtain quality data from these raw signals, 
the two output pulses from the laser have been interleaved 
to prevent crosstalk and potential electronic noise, which can 
contaminate the first few range bins of the retrieval. Inter- 
leaving the pulses also puts an offset into the system that 
must be accounted for. This offset is important because of 
the sensitivity of the DIAL equation (5) in the retrieval near 
the surface. This offset was empirically determined using a 
photodiode and an oscilloscope, but is also shown to be cor- 
rect in the final retrieved results. 

Due to the dynamic nature within the PBL and lower free 
troposphere, it may be necessary to resolve fine scale mixing 
features or a layer of ozone aloft at an upper altitude. De- 
pending on the research application, the temporal resolution 
is selected and the raw data sets are averaged and is defined 
by the number of data sets that are averaged. For example, 
a 10 minute temporal resolution corresponds to the running 
average of 30 raw signal datasets. The time averaged data is 
then background subtracted using altitudes where the slope 
of the background is close to zero, in order to assume a lin- 
ear background correction. This region is typically greater 
than 30 km. 

After these corrections, it may be necessary to manually 
remove data sets if they are believed to be contaminated by 
clouds, which is typically quantified by a large rise in lu- 
minosity counts at specific altitudes in the individual detec- 
tors. If the cloud optical thickness is fairly low and multiple 
data sets are averaged together, information may be available 
above the cloud layer. 

The retrieval algorithm must also correct for the nonpar- 
alyzable dead time correction of the photomultiplier tubes 
(Keckhut et al. (2004)), which has been empirically deter- 
mined to be between 4 and 7 ?rs for most lidar return sig- 
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nals. These values are quantified by comparing the lidar re- 
turn signal to a model atmosphere from a standard model or 
co-located sonde data in the Low channels. Subsequent lidar 
returns from the Middle and High channels are then corrected 
to match the overlapping regions of the Low signal returns. 730 
This correction can be written as, 


C t = 


C„ 


1 — C m T d ’ 


(9) 


where the true photon count rate ( Ct ) can be expressed as a 
function of the measured count rates (Cm) and a dead time 
(Td) parameter (Lampton and Bixler (1985)). This correction 
is critical in the Low channels where the PMT signals cover a 

735 

large dynamic range and the derivative in the DIAL equation 
is large. Signal Induced Bias (SIB) is also quantified, but the 
channel merge altitudes are selected in a way such that SIB 
is not a dominant correction. 

Referring to equation (5), the retrieval relies on the deriva- 
tive of the natural logarithm of the ratio of backscattered 
laser powers. The finite impulse response (FIR) Savitsky- 
Golay (SG) differentiation filter has been used to produce 
the required first order derivative. The SG filter is a gener- 
alized running average with filter coefficients determined by 
an unweighted linear least-squares regression and a 2nd or- 

740 

der polynomial model applied to the derivative. The second 
order was chosen instead of a third or fourth order because 
it was less likely to pick up extreme noise in the derivative. 
The final vertical resolution of the data is determined using 
the full width at half maximum (FWHM) of the steady state 

745 

SG filter coefficients associated with the filter window size. 

Due to a known temperature dependence of the ozone ab- 
sorption cross sections (Figure 1, Malicet et al. (1995)), it 
is necessary to get an accurate atmospheric temperature pro- 
file, either from a co-located sonde launch of from a stan- 
dard model atmosphere. Because the ozone absorption tem- 
perature dependence is not known continuously, a 2nd order 
polynomial interpolation is utilized. The Rayleigh molecu- 
lar extinction coefficients are also corrected for atmospheric 
number density values provided or derived from a co-located 
sonde launch or from a standard model. After these cor- 
rections are made, the first estimate of the ozone profile, 
No 3 (r)i, is computed. 755 

It is necessary to correct the computed ozone profile in 
instances of strong aerosol influence, which are particularly 
prevalent within the PBL. The aerosol extinction term in 
equation (7) depends not only on the wavelength of the in- 
teracting photons, but also on the aerosol size and type. To 
quantify this, an iterative algorithm has been used to deter- 
mine the aerosol backscatter coefficient, which when cou- 
pled with an assumed lidar ratio, will yield an estimate of 760 
the aerosol extinction coefficient. The lidar ratio, S, is the 
extinction to backscatter ratio. 


S = a(r)/P(r). 


(10) 


This procedure, which utilizes the offline, less absorbing 
wavelength of 299 nm, begins with an initial estimate of 
the aerosol volume backscatter at a reference height ( h re f ) 
and constant lidar ratio which are adapted from Kuang et al. 
(2011) and Browell et al. (1985). 

Rearranging the final DIAL equation (5) with the aerosol 
correction term, it is possible to solve for the first estimate of 
the volume backscatter coefficient due to aerosol as, 


Paerij") 1 


exp (A + B) X 


(Prnol (t') T P aer (^re /)) 
(r + Ar ) 2 


Pmoi{r) 

( 11 ) 


where 


JWr + A,-) 
-* 299 v ) 


and 


( 12 ) 


Bi = -2No 3 (r)i<Jo 3 Ar - 2Ar(a moi (r) - Sp ae r{h re} )). 

(13) 

The previously obtained non-aerosol corrected ozone 
number density, No 3 (r) 1 , is used for the first guess of 
the ozone profile in order to calculate B\. The remain- 
ing molecular components are calculated from a standard 
model atmosphere or from a recent sonde launch. The aerosol 
backscatter coefficient at the reference height that is used is 
/3aer{h re f) = 1.67 x 10 ~ 4 m~ 1 sr~ 1 , h re f is 12 km and the 
lidar ratio is assumed to be S = 60sr -1 . At this point, equa- 
tion (11) has been solved for the first complete profile of the 
aerosol backscatter coefficient, fl aer (r) i. 

For the second iteration, a more accurate value of the 
aerosol backscatter coefficient can be computed by making 
the substitution of 


Sfiaerihref ) = S(P aer (r)i + Paev{r + Ar) 1 )/2. (14) 

This can be solved for the new aerosol backscatter coeffi- 
cient, Paer(r) 2 - 

This iterative process continues, making the necessary 
substitution of 

P aer 0)i+l + Paer + Ar ) i+ i = p aer (r)j + P aer (r + A?’); , 

(15) 

until a stable solution is reached, where the difference in iter- 
ations of the profile is negligible. The final aerosol backscat- 
ter term is then inserted into the full DIAL equation (6) and 
(7) for the final corrected ozone mixing ratio profile. In the 
test case below, the calculated correction in the final reported 
ozone mixing ratios were <3 ppbv in the Low channels and 
<1 ppbv in the Middle and High channels. 
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The final aerosol corrected ozone profile is then joined to 
form one continuous hybrid profile. The profile join regions 
are chosen to maximize SNR, which is at its highest in the 
lower range of each profile and significantly deteriorates ats 2 o 
the top range of each profile. With the current configuration 
for a 10 minute temporal running average and a fixed SG 
filter window size, the vertical resolution of the Low channels 
is 275 meters, the Middle channels is 525 meters and the 
High channels is 750 meters. 825 

5 Comparison with Beltsville, MD Ozonesonde [25 Oc- 
tober 2013] 

830 

As mentioned before, the Mid-Atlantic region displays a 
moderate amount of seasonal ozone variability with the high- 
est concentrations occurring in the late July time period 
(Lehman et al. (2004)). Because of this, during the fall and 
winter seasons, the TROPOZ focuses its observations on 835 
monitoring ozone in the free troposphere and the aforemen- 
tioned source contributions, such as STE. While taking mea- 
surements on 25 October 2013, an unusually low tropopause 
was observed, resulting in the exchange of high amounts of 
ozone between the stratosphere and the lower free tropo- 84o 
sphere. An ozonesonde was launched that day by the Howard 
University Beltsville Center for Climate Systems Observa- 
tion, which serves as the validation source for the TROPOZ 
ozone profiles. The launch site (39.05°N, 76.88°W) is ap- 
proximately 8 km from the lidar site which is close enough 845 
to assume similar, but not identical tropospheric micromete- 
orology in the dynamic daytime PBL. 

A 10-minute retrieved lidar ozone profile (17:44:1754 
UTC) together with the ozonesonde profile (17:44 UTC 
launch) are shown in Figure 5. The retrieved range of the sso 
lidar profile is from 400 m to 12 km AGL. Both the lidar and 
the sonde show a PBL height near 1.5 km which is based 
on the slight ozone mixing ratio gradient, the strong relative 
humidity gradient, and the temperature inversion. Above the 
PBL, both profiles have values in the range 30-40 ppbv un- 855 
til the strong ozone gradient near 6.5 km. Both the sonde and 
the lidar capture this strong ozone reservoir, where the differ- 
ence in the altitude of the reservoir may be attributed to the 
distance between sites based on the fact that the ozonesonde 
profile still resides within the uncertainty of the retrieved pro- 860 
file. This ozone reservoir, which was around 150-200 ppbv, 
had a higher ozone mixing ratio than most summer days and 
previous test cases did not reach these values until 12 km, 
which emphasizes the regional and seasonal uniqueness of 
this test case. The uncertainty bars shown in the retrieved 865 
TROPOZ ozone profile are dominated by the statistical un- 
certainty of photon counting and SNR, which are discussed 
more in Section 6. 

The TROPOZ profile differs from the sonde profile above 
8 km; an explanation for this can be seen in Figure 6. This im- 870 
age shows the GPS trajectory of the sonde for its first 10 km 


of ascent overlain on a map of the area. The color of the tra- 
jectory track represents the altitude of the sonde. The sonde 
was launched at Beltsville, MD and was closest to GSFC in 
the ascent at an altitude near 3 km, which also implies that 
the difference in ozone mixing ratio in the first 2 km may 
be due to the distance between sites. As the sonde continued 
along this trajectory, it entered a marine environment. The 
eastern end of the trajectory, where the sonde had achieved 
10 km altitude, is near Edgewater, MD, nearly 30 km dis- 
tant from the lidar site. Because the lidar is in an Eulerian 
frame of reference in the inland urban environment, it was 
not actually resolving the same air mass that is present at the 
beginning or end of the Lagrangian frame of reference of the 
sonde. This highlights one of the inherent complications in 
using ozonesondes as a validation source for lidar. Because 
the lidar is stationary, the resolved ozone profiles are more 
representative of the ozone profile at the observation site. 

The ozone mixing ratio, relative humidity and tempera- 
ture are plotted from the ECC and radiosonde in Figure 7. 
The ozone mixing ratio is fairly constant in the first 6 km 
with low variability until the large ozone gradient at 6.5 km. 
Above that point there are large oscillations that imply a dy- 
namic system of ozone in this region. The relative humid- 
ity profile shows a strong gradient at the top of the PBL but 
above that the air slowly gets drier until about 6.5 km where it 
also undergoes a strong gradient. The relative humidity also 
shows some small oscillations that correspond to the same 
altitudes as the dynamic ozone oscillations. The temperature 
profile exhibits a linear decrease up to 6.5 km above which 
the temperature remains nearly constant for 2 km. This, to- 
gether with the observation of a dry air mass above 6.5 km, 
indicate an unusually low tropopause height near 6.5 km. 
This implies that this air mass is most likely of stratospheric 
origin and can be considered an ozone intrusion event and 
can be classified as STE (Reiter (1975)). 

A 500-millibar (mb) height contour national map for 1 1 :00 
UTC 25 October 2013, courtesy NOAA Daily Weather Map, 
for the meteorological events that surrounded the test case is 
shown in Figure 8. The chart shows height contours (solid 
lines), temperatures (dashed lines) and winds (arrows). The 
height contours show the height of the 500-mb pressure level 
in dekameters above sea level, the isotherms, or lines of con- 
stant temperature, are shown in degrees Celsius, and the ar- 
rows show the wind direction and speed at the 500-mb level. 
The map shows a persistent low pressure system occurring in 
the Mid-Atlantic and Ohio River Valley regions just under- 
neath the tropopause height around 6.5 km. This disturbance 
is exemplified by the deviation from the trough of closely 
spaced geopotential height contours (solid black lines) and in 
the temperature contours (dashed black lines). Also, notice 
the strong increase in northwesterly winds associated with 
the steep pressure gradients that have occurred just below the 
tropopause around 6.5 km. This wind pattern is most likely 
the Polar jet stream (Viezee et al. (1983)), which is typically 
stronger in winter and fall months and was the reason why 
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the ozonesonde in Figure 6 moved so quickly in the zonal 
direction away from the launch site as it got nearer to the 
tropopause height. The low pressure system was present in 
national maps as early as 23:00 UTC 24 October 2013 and 930 
slowly subsided and evolved throughout 25 October 2013 as 
can be seen in the lidar curtain plots in Figure 9a and 9b. 

Figure 9a shows the time series from 15:19-22:18 UTC 
(11:19 - 18:18 EST) with the overlaid ozonesonde data corre- 
sponding to the red arrow point at 17:44 UTC. This image is 935 
created from a 10 minute running average of the data with all 
altitude channels merged together to form a complete profile 
of the atmosphere from 400 meters to 12 km. It is possible to 
see the strong gradient in ozone near 7 km and large ozone 
concentrations that dissipate mainly through titration during 940 
the afternoon and early evening. The upper level cold and 
dry air pockets behind the persistent low pressure system in- 
ject a large ozone reservoir (150-200 ppbv) into the free tro- 
posphere, effectively lowering and maintaing a tropopause 
height near 6.5 km. The in situ surface ozone mixing ra-945 
tios from the GSFC TROPOZ DIAL trailer (Thermo Scien- 
tific 49i, 1 minute resolution) are reported in the first bin in 
the time series below the retrievable altitude. These values 
range from 18-24 ppbv at the surface for the entire experi- 
ment, which implies at ground level the ozone was behav-950 
ing like a typical early winter Mid-Atlantic regional atmo- 
sphere. Based on the northwest direction and strong wind 
speeds shown in Figure 8 , there was dilution and ventilation 
of ozone precursors which resulted in a low amount of of 
ozone near the surface from local pollution. 

Around 19:00 and 21:00 UTC, thin clouds drifted through 
the field of view during the observation, the effects of which 

^ 955 

can be seen as white space in the figure. Because the 10 
minute observation is a running average over 30 profiles, and 
because the clouds were thin and did not form broad cover, 
retrieval information can be seen above the clouds in the cur- 
tain plot. 

The curtain plot in Figure 9b shows the same data set with 
a temporal resolution of 1 minute from 4 to 10 km. This is 
a high enough resolution to begin to pick out some of the 
small scale features in the ozone curtain plot. With an run- 960 
ning average of just 3 profiles, retrieval above the clouds is 
lost, hence, the vertical white lines that appear in Figure 9b. 
Averaging less profiles does increase the uncertainty in the 
measurement and the retrieved ozone profile at 1 minute res- 
olution can only reach 10 km (as opposed to 12 km with a 10 965 
minute resolution) before the SNR is too small yield useful 
information. 

This plot focuses attention on the upper atmospheric con- 
ditions in order to represent the ability to monitor the chang- 
ing dynamics of the ozone reservoir on a very short time 970 
scale. It is possible to see small pockets of clean air which 
implies that there was some dynamic mixing of clean tropo- 
spheric air and high ozone stratospheric air. There is also an 
interesting wavelike feature in the ozone concentration to- 
wards the end of the run (20:30-22:18 at 6.5 km) that may 975 


have been the cause for the ozone reservoir to split and be- 
gin to both subside towards the surface and to rise back to 
the upper free troposphere. The ozone free region between 
two ozone enhanced layers is certainly characteristic of a 
tropopause fold event (Viezee et al. (1983)), but without the 
meteorological parameters from a radiosonde, it is difficult to 
confirm the fold. Dilution of this ozone may have occurred in 
the atmosphere during the times at which the TROPOZ was 
not taking measurements, especially with high winds push- 
ing the air mass out towards the Chesapeake Bay. 

Surface monitors did show an increase in ozone on 26 Oc- 
tober 2013; but, with the passing front bringing warm, sunny 
weather, it is challenging to determine whether the surface 
ozone source was due to the dynamic STE or created by pho- 
tochemistry within the PBL. 

Due to the unusually high amount of ozone in the atmo- 
sphere during this test case, the usable return signal was ab- 
sorbed much more strongly than a typical day. Because of 
this, ozone mixing ratios are only reported to 12 km. For 
conditions of nighttime observations and longer temporal av- 
erages, ozone retrievals can be attained at a higher altitude. 
Although this plot emphasizes the upper atmospheric condi- 
tions for this regionally important test case, it provides evi- 
dence that the GSFC TROPOZ DIAL will be able to provide 
high resolution ozone observations within the dynamic PBL. 

6 Accuracy of the DIAL measurement 

The GSFC TROPOZ DIAL system detects individual pho- 
tons through the use of photomultiplier tubes. The signal col- 
lected by these PMTs follows Poisson statistics (Megie et al. 
(1985), Papayannis et al. (1990)) and the statistical uncer- 
tainty of the ozone number density can be calculated as 


ejv „ 3 (r) 


1 I Ps,\ + Pb, A + Pd,\ 

2N 03 AaoAr\j ^ 7 “ 


(16) 


where P s .\, Pb,x, and P^.a are the atmospheric backscattered 
signal, background radiation, and dark counts of the detector 
at wavelength A. The vertical resolution, which is based on 
the SG filter window size, is denoted as Ar and the differen- 
tial ozone absorption cross section is denoted as AtJo 3 - The 
statistical uncertainty is related to the square root of the to- 
tal PMT counts, both those that are relevant to the retrieval 
of ozone number density and those that are counts due to 
systematic uncertainties. By integrating profiles for a longer 
duration, the backscattered signal term P s \ becomes much 
larger than the P;, a and Pd, a terms. For this reason, the tem- 
poral resolution is inherently built into the statistical uncer- 
tainty of the system and averaging many data sets is benefi- 
cial to the resultant uncertainty in the system. 

By increasing the vertical resolution of the profile, which 
is in the denominator of equation 16, the statistical uncer- 
tainty in the measurement decreases. Because the vertical 
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resolution changes with the retrieval range and has different 
values for different channel pairs, the resultant uncertainty 
profile exhibits changes at the altitudes at which the chan-1030 
nels are merged. This can be seen in the uncertainty bars in 
Figure 5. 

Aside from the statistical uncertainty, there are uncertain- 
ties associated with the spectral dependence of other inter- 
fering species, such as aerosols, non-ozone interfering gases , 1035 
and the Rayleigh atmosphere. As mentioned in the Data Pro- 
cessing portion of this paper, an iterative aerosol correction 
was applied to this test case, which resulted in a <3 ppbv 
(< 12%) ozone correction in the Low channels and a < 1 ppbv 
(<4%) correction in the Middle and High channels. 1040 

Further uncertainties arise from interfering gases that have 
similar absorption cross sections at the DIAL wavelengths 
(Keckhut et al. (2004)). An uncertainty of 1.5% in overall 
retrieved ozone has been calculated due to absorption at the 
DIAL wavelengths by the oxygen dimer (02-02) interaction 045 
(Fally et al. (2000)), a -1.0% uncertainty due to sulfur diox- 
ide (S02) (Stark et al. (1999)), and a -0.5% uncertainty due to 
nitrogen dioxide (AY^KBogumil et al. (2003)). The uncer- 
tainty in the Rayleigh extinction correction is typically <1% 
due to the dependence on the atmospheric number density,io5o 
which may be calculated from a standard model atmosphere 
or recent radiosonde. Similarly, the temperature dependence 
of the ozone absorption cross sections should contribute less 
than 3% to the uncertainty in the final ozone number density 
(Leblanc et al. (1998)). 1055 

Table 3 lists the various uncertainties due to individual 
photon counting, interfering gases, the temperature depen- 
dence of the ozone cross section, and the saturation correc- 
tion for the test case discussed in the previous section. Other 
than the uncertainty due to the presence of aerosols, and toio6o 
a lesser extent those due to Rayleigh extinction and satura- 
tion correction, these values can be taken as representative of 
the uncertainties in ozone retrieval for the GSFC TROPOZ 
DIAL instrument. The uncertainties are largest in the Low 
channel due to the presence of aerosols and the complete usao65 
of the saturation correction. In the Middle and High channels, 
there are fewer aerosols present and there is enough over- 
lap between the Low and Middle and between the Middle 
and High channels to join them at altitudes such that regions 
where saturation is non-negligible can be excluded from thao7o 
overall profile. The reported resultant uncertainty, added in 
quadrature, from 0-1.5 km is 16-19%, 1.5 -3 km is 10- 18%, 
and from 3 km to 12 km is 1 1 - 25%. 
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The GSFC TROPOZ DIAL has been utilized to detect tro- 
pospheric ozone, a very important greenhouse gas for cli- 
mate and health studies. Due to complex tropospheric chem- 
istry and various free tropospheric ozone sources, modeling 
these lower levels of the atmosphere can be extremely difioso 


ficult and empirical observations of ozone may be far more 
useful. The necessary theory and initial optical layout of the 
GSFC TROPOZ DIAL in order to retrieve and better char- 
acterize the vertical distribution of tropospheric ozone in the 
Baltimore-Washington D.C. area has been established. A re- 
gionally important validation test case of an unusually low 
tropopause, resulting in the exchange of high amounts of 
ozone between the stratosphere and the lower free tropo- 
sphere, or STE, is presented and analyzed. 

Measurements of the STE are presented with a temporal 
resolution of 10 minutes and 1 minute with a maximum alti- 
tude of 12 km. Because of the optically thick stratospheric 
ozone layer, a ground based detection method, which can 
detect tropospheric ozone on a continuous, high resolution 
scale is advantageous. The GSFC TROPOZ DIAL system 
has shown to be a less expensive and more robust alternative 
to previous lidar systems with the implementation of Raman 
cells, which have had no observable drop in pressure since 
they were initially filled in lune 2013. The data processing 
algorithm has been introduced and has shown to have agree- 
ment with a co-located ozonesonde measurement. The re- 
ported resultant uncertainty from 0-1.5 km is 16-19%, 1.5 
-3 km is 10 - 18%, and from 3 km to 12 km is 1 1 - 25%. 

Future plans for the GSFC TROPOZ DIAL are to obtain 
ozone profiles as close to the surface as possible (< 100 m 
AGL), which may be obtainable with ozonesonde launches 
directly at the lidar site. As mentioned previously, the Raman 
cells were filled to a pressure which left a residual amount 
of 266 nm for atmospheric transmission. Because this wave- 
length is strongly absorbed by ozone it may only be able 
to yield useful information from the surface to the top of 
the PBL, but may have enough signal to characterize and 
probe the dynamic PBL. The advantage of adding 266 nm 
to this system would be to get a more robust aerosol correc- 
tion method using a dual-DIAL wavelength retrieval. 

Depending on the atmospheric event the system is in- 
tended to monitor, averaging of the data may be minimized 
to detect natural variability and small ozone fluctuations or 
maximized to reach higher altitudes for climatology pur- 
poses. The GSFC TROPOZ DIAL is a new and unique tool 
for the air quality community in the Baltimore-Washington 
D.C. region because it can provide a continuous, high reso- 
lution and novel perspective of ozone pollution episodes. For 
these reasons, the GSFC TROPOZ DIAL is a very valuable 
instrument for quantifying the vertical distribution of ozone, 
especially for distinguishing layers of ozone aloft and STE 
events that surface monitors do not have the ability to detect. 
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Table 1: Vibrational Raman frequencies and first Stokes wavelengths obtained from the stimulated Raman scattering of 
Nd:YAG fourth harmonic in the Raman active gases. [Source: Measures (1983)] 


Select Gas Vibr. Raman Frequency 1st Stokes 

Hydrogen 4160.0 cm -1 299.1 nm 

Deuterium 2986.1 cm -1 288.9 nm 
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Table 2: Hardware Description for the GSFC TROPOZ DIAL system for the transmitter, receiver components and the data 
aquisiton system. 


System Component 

299 (“off’) Channels 

289 (“on”) Channels 

Wavelength Source 

Nd:YAG 

Nd:YAG 

Pulse Width 

8 ns 

8 ns 

Beam Divergence 

0.42 mrad 

0.44 mrad 

Repetition Rate 

50 hz 

50 hz 

Pump Power 

32 mJ/pulse (1.6 W) 

42 mJ/pulse (2.1 W) 

Raman Cell Partial Pressures [psig] 

200 H 2 /600 He 

300 D 2 /400 He 

Gas Purities 

99.9999 % H 2 

99.999 % D 2 

Raman Cell Efficiencies 

53% 

26% 

Transmitted Raman Pulse Energy 

17 mJ/pulse (0.85 W) 

1 1 mJ/pulse (0.55 W) 

Interference Filters 

Materion Precision 

Materion Precision 

Center Wavelengths 

299.1 ± 0.25 nm 

288.9 ± 0.25 nm 

FWHM 

<1.2 nm 

<1.2 nm 

Low Channel [0-1.5 km] 

Telescope Diameter [cm] 

2.5 cm 

2.5 cm 

FOV [mrad] 

4.5 

4.5 

Detector (Hamamatsu) 

7400P-03 

7400P-03 

Gates [km] 

0.15 

0.15 

Neutral Density Transmittance [%] 

5 

10 

Calculated Full Overlap Altitude [km] 

0.105 

0.105 

Middle Channel [1.5 -5 km] 

Telescope Diameter [cm] 

2.5 cm 

2.5 cm 

FOV [mrad] 

4.5 

4.5 

Detector (Hamamatsu) 

7400P-03 

7400P-03 

Gates [km] 

0.415 

0.415 

Neutral Density Transmittance [%] 

25 

50 

Calculated Full Overlap Altitude [km] 

0.105 

0.105 

High Channel [3-12 km] 

Telescope Diameter [cm] 

45 cm 

45 cm 

FOV [mrad] 

1.0 

1.0 

Detector (Hamamatsu) 

7400P-03 

7400P-03 

Gates [km] 

3.015 

3.015 

Neutral Density Transmittance [%] 

100 

50 

Calculated Full Overlap Altitude [km] 

0.462 

0.462 


Data Acquisition System Licel Transient Recorder (Model TR20-80) 

300 Mhz Maximum 
Photon Counting Rate 
16-bit, 50-ns range resolution 
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Table 3: A list of the quantified uncertainties associated with the GSFC TROPOZ DIAL instrument during the test case on 25 
October 2013 


Low Channel [0-1.5 km] 

Middle Channel [1.5-5 km] 

High Channel [3-12 km] 

Statistical Uncertainty 

1-10% 

1-16% 

1-24% 

Aerosol Correction 

< 12% 

<2% 

<2% 

Rayleigh Correction 

< 1% 

< 1% 

< 1% 

S 0 2 

-1% 

-1% 

-1% 

NO2 

-0.5% 

-0.5% 

-0.5% 

O2 - Dimer 

1.5% 

1.5% 

1.5% 

Temperature Dependence of A<to 3 

3 % 

3% 

3% 

Saturation Correction 

10% 

5% 

3% 

Final Reported 

16-19 % 

10-18% 

11-25% 
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Ozone Absorption Cross Section 



Wavelength (nm) 


Fig. 1: The spectral dependence of the ozone absorption cross section in the UV. The more absorbing wavelength, at 288.9 nm, 
is called the “on” wavelength and the less absorbing wavelength, at 299.1 nm, is called the “off' wavelength. [Source: Malicet 
et al. (1995)] 
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Photon Conversion Efficiency [Hydrogen & Helium] 



Fig. 2a: Photon conversion efficiency of hydrogen as function of total Raman cell pressure. The experiment was done with 
multiple partial pressures of hydrogen and helium, in which the the most efficient pressure combination of S 1 was chosen for 
the GSFC TROPOZ DIAL. The pressures indicated in the legend are for the partial pressures of the Raman active gas. The 
uncertainty bars represent one standard deviation in the measurements. 


Photon Conversion Efficiency [Deuterium & Helium] 



Total Cell Pressure (psi) 


Fig. 2b: The experiment repeated for the multiple pressure combinations in the deuterium Raman cell. 
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Fig. 3: The experimental schematic that includes the 50 hz Nd:YAG laser, steering optics, Raman cells, near field smaller 
telescopes, main telescope and detection packages. 
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Fig. 4: The 13 m trailer that houses the GSFC TROPOZ DIAL system. The hatch doors are open to transmit the DIAL 
wavelengths into the atmosphere. 
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GSFC TROPOZ DIAL [10 Min Avg.] 
25-Oct-2013 17:44 (UTC) 



Fig. 5: The profile comparison for the GSFC TROPOZ DIAL retrieval (10 minute resolution) and the ozonesonde launch 
at 17:44 UTC on 25 October 2013 at Beltsville, MD. [Source: Howard University Beltsville Center for Climate Systems 
Observation] 
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Fig. 6: Ozonesonde GPS coordinates during the ozonesonde flight. The color represents the altitude of the ozonesonde with 
the first red dot corresponding to Beltsville, MD. The TROPOZ site at GSFC is near the green section of the trajectory. When 
the sonde is above 7 km (light and dark blue) it is 16 km away from the lidar and at a altitude of 10 km (magenta) it is near 30 
km away from the lidar. The altitude legend is on the right side of the figure. [Source: Howard University Beltsville Center for 
Climate Systems Observation; Map Source: Google Maps] 
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Fig. 7: The radiosonde and ECC in situ data for ozone mixing ratio, relative humidity, and temperature plotted for the test case 
on 25 October 2013. [Source: Howard University Beltsville Center for Climate Systems Observation] 
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Fig. 8: A 500-millibar (mb) height contour national map for 11:00 UTC 25 October 2013, courtesy NOAA Daily Weather 
Map, for the meteorological events that surrounded the test case is shown. The figure shows height contours (solid lines), tem- 
peratures (dashed lines) and winds (arrows). The height contours show the height of the 500-mb pressure level in dekameters 
above sea level, the isotherms, or lines of constant temperature, are shown in degrees Celsius, and the arrows show the wind 
direction and speed at the 500-mb level. [Source: NOAA Daily Weather Map] 
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GSFC TROPOZ DIAL 
25-Oct-2013 15:19 - 22:18 UTC 
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Fig. 9a: Time series plot of the ozone mixing ratio [ppbv] from the GSFC TROPOZ DIAL between 1 5 : 19-22: 1 8 UTC (11:19 - 
18:18 EST) at a 10 minute temporal resolution with the co-located ozonesonde launch corresponding to the red arrow point at 
17:44 UTC. The surface ozone monitor mixing ratios (1 minute resolution) are plotted in the lowest bin of the figure. 
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Fig. 9b: Time series plot of the ozone mixing ratio [ppbv] from the GSFC TROPOZ DIAL between 15:19-22:18 UTC (11:19 
- 18:18 EST) at a 1 minute temporal resolution with the co-located ozonesonde launch corresponding to the red arrow point 
at 17:44 UTC. This plot focuses attention on the upper atmospheric conditions in order to represent the ability to monitor the 
changing dynamics of the ozone reservoir on a very short time scale. 
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